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Abstract 

Posttranslational modification of proteins is key in transmission of signals in cells. 
Many signaling pathways contain several layers of modification cycles that mediate 
and change the signal through the pathway. Here, we study a simple signaling 
cascade consisting of n layers of modification cycles, such that the modified protein 
of one layer acts as modifier in the next layer. Assuming mass-action kinetics and 
taking the formation of intermediate complexes into account, we show that the steady 
states are solutions to a polynomial in one variable, and in fact that there is exactly 
one steady state for any given total amounts of substrates and enzymes. 

We demonstrate that many steady state concentrations are related through ratio- 
nal functions, which can be found recursively. For example, stimulus-response curves 
arise as inverse functions to explicit rational functions. We show that the stimulus- 
response curves of the modified substrates are shifted to the left as we move down 
the cascade. Further, our approach allows us to study enzyme competition, seques- 
tration and how the steady state changes in response to changes in the total amount 
of substrates. 

Our approach is essentially algebraic and follows recent trends in the study of 
posttranslational modification systems. 

Keywords: Mass action kinetics, Sequestration, Futile cycle, Post-translational 
modification, Stimulus-response, Rational function 

1 Introduction 

Posttranslational modification of proteins is one of the principal mechanisms by which 
signals are transmitted in living cells. In particular, covalent modification by phosphory- 
lation is wide-spread and provides one of the basic modes by which a signal is mediated 
through a pathway. Classical signaling pathways typically contain a cascade of post- 
translational modification cycles (also called futile cycles) where the activated protein 
in one layer acts as the modifier enzyme in the next layer. In order to understand the 
biological importance of this sophisticated signaling mechanism and to predict the be- 
havior of the system, a number of theoretical studies have focused on determining the 
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system's dynamics and steady states. In the most discussed case, the MAPK cascade, 
the presence of bistability and oscillatory behavior has been revealed for different choices 
of system parameters [SI E3] • 

The basic features of general signaling cascades can be elucidated from the study 
of signaling cascades with n layers and one cycle of post-translational modification at 
each layer. These cascades are part of many pathways (examples include the cAMP 
pathway [6, page 342], quorum sensing pathways [33] and the intrinsic cascade in the 
blood coagulation pathway [19]) and have been extensively studied from different math- 
ematical points of view: n = 1, e.g. [31 [26], n = 1,2, e.g. [9[ [TO], [18], and arbitrary n, 
e.g. [U HU EU EI] . Here, we use the classical model of Michaelis and Menten in which 
the enzyme-substrate complex is formed reversibly and the dissociation of the complex 
into a modified substrate and a free enzyme is assumed to be irreversible. Further, the 
modifier, which in the case of phosphorylation typically would be ATP, is assumed to 
be in constant concentration and hence not included explicitly in the model. In this we 
follow much previous work, e.g. [9] fTT| IHTj . 

In order to simplify the mathematics of the Michaelis-Menten model, additional re- 
quirements are often imposed; e.g. the formation of intermediate complexes is ignored 
under the quasi steady state assumption [H [TH [23] or the total amount of substrate is 
assumed to be much higher than those of the (de)modification enzymes [31] • In the first 
case, this has the consequence that some biological mechanisms, like sequestration, might 
be overlooked, while the second case is not, in general, considered biologically valid [3J. 
Additionally, conclusions are often drawn from parameter sampling and are not based 
on a rigorous mathematical analysis of the model [23] . 

We focus on understanding the steady state behavior of a signaling cascade with n 
layers and one cycle of post-translational modification at each layer. Using mass-action 
kinetics and without imposing further assumptions on the Michaelis-Menten model, we 
show that the steady states arise as the solutions to one polynomial in one variable. This 
allows us to conclude that these systems have exactly one steady state for fixed total 
amounts of substrates and enzymes. This is achieved using an induction argument that 
reflects the modularity of the cascades. 

Our approach makes it possible to study aspects of the system in detail without rely- 
ing on simulation. For example, the stimulus-response curve arises as the inverse function 
of an explicit rational function (i.e. a quotient of polynomials in one variable). Several 
other rational functions are found that describe the species concentrations at steady 
state; these functions provide a clear picture of how the steady state re-accommodates 
to variations in initial concentrations and provide insight into sequestration. In addition, 
our approach enables us to explore the parameter space in more detail than previously 
and to divide the space into regions that exhibit different steady state behaviors. 

The simplicity exposed by the algebraic approach is the key component in this work 
and follows the line of recent works [HI [29J [30] , as well as some older approaches [2] . All 
rational functions can be explicitly computed. In addition, their inverses can be obtained 
using any program that allows symbolic algebraic computations like Mathematica™ . 

The article is divided into two main sections. In the first the mathematical theory 
is developed. We show that the steady states can be described as the solutions to a 
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polynomial in one variable only and that there is only one steady state for any initial 
total amounts of substrates and enzymes. Subsequently, we derive rational functions 
that relate the concentrations of chemical species at steady state. In the second section, 
we use these results to discuss various topics of biological interest such as enzyme com- 
petition, stimulus-response curves, and how steady state concentrations vary in response 
to changes in the total amount of substrates. We end with some concluding remarks. 
Proofs are given in the appendix. 

Acknowledgements. This work was done while EF was visiting Aarhus University in Spring 
2010. It was completed while EF and C W were visiting Leipzig University in July 2010. EF wishes 
to thank the members of the Bioinformatics Research Centre in Aarhus for their hospitality. EF 
is supported by the Spanish Ministerio de Ciencia e Innovation by project MTM2009-14163-C02- 

01 and a travel grant from the Universitat de Barcelona. MK, LNA and CW are supported by 
the Danish Research Council and CW by the Humboldt Foundation. All authors are supported 
by a travel grant from the Carlsberg Foundation, Denmark. Joaquim Puig is thanked for useful 
discussion. 

2 Mathematical results and characterization 
2.1 Preliminaries 

By a rational function f(x) over a field K we understand a function that is given as a 
quotient of polynomials f(x) = p(x)/q(x) with coefficients in K. Assuming that it is 
in reduced form, that is, there are no common zeros between p{x) and q(x), then the 
only discontinuities of the function are given by the zeros of the denominator. Note 
that two reduced forms p±(x) / q±(x) and P2{x)/q2{x) of the same rational function satisfy 
Pi(x) = \p2^x) and q\{x) = Aq^x) for some A E K. If not otherwise stated, we consider 
rational functions over R. 

If S is a finite set, M[S] denotes the ring of polynomials in S and M(S) denotes the 
extension field of M. consisting of quotients of polynomials in S, that is, the fraction field 
of M.[S] (see e.g. p2])- There is an evaluation map from M[S] to M. obtained by assigning 
a real value to every element of S. 

Let IR+ denote the set of positive real numbers (excluding zero) and let M + denote 
the set of non- negative real numbers (including zero). An increasing (resp. decreasing) 
continuous function f{x) defined on an interval / C ]R is a continuous function such that 
f(x) > f(y) (resp. f(x) < f(y)) for x,y £ I satisfying x > y, i.e. we consider strictly 
increasing/decreasing functions. 

We will make several uses of the following result: 

Lemma 2.1. Let f(x) = p{x)/q{x) be a rational function of the variable x, that is, a 
quotient of real polynomials p(x) and q{x) in x. Let L CM be an open interval. If q(x) 
has no zeros in I, then f(x) is a continuous function in I. In addition: 

(i) If f'(x) > for all x £ I, then f(x) is an increasing function in I. 
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(ii) If f'{x) < for all x £ I, then f(x) is a decreasing function in I. 

If either (i) or (ii) are satisfied, then by the Inverse Function Theorem, f admits a con- 
tinuous inverse function g: im(/) — > I. The function g is increasing (resp. decreasing) 
if f is. 

If q(x) is a real polynomial of degree to, then q(x) has at most to real zeros (it has 
exactly to complex zeros when counted with multiplicity). If x\ < ■■■ < x r are the 
different real zeros of q(x), ordered increasingly, then f(x) is a continuous function in 
each of the r + 1 intervals (— oo, xi), (x r , +oo) and (xj, Xj+i) for i = 1, . . . , r — 1. 

2.2 One-site PTM signaling cascades 

We consider signaling cascades with n layers and one posttranslational modification 
(PTM) cycle at each layer, as illustrated in Figure [TJ The chemical species involved in 
each cycle are the unmodified substrate Sf, the modified substrate Sf, the demodification 
enzyme F{, the modification enzyme Ei = S}_ 1 for i = 2, . . . ,n, and the intermediate 
enzyme-substrate complexes Y® and Y^. That is, in each layer the modification enzyme 
is the activated substrate of the previous layer. The enzyme of the first layer, E\, is not 
a substrate in any other layer. For convenience we put Sq = E\ = E. In the following, 
signaling cascades of this kind are called one-site PTM cascades, or just cascades for 
short. 

Although we consider PTM cycles of any type, we use the nomenclature of modifica- 
tion by phosphorylation and call the catalyzing enzymes for kinase (modification) and 
phosphatase (demodification). 

The system is specified by the set of chemical reactions shown in Figure [TJ The 
enzyme mechanism is assumed to follow the classical model of Michaelis and Menten, in 
which an enzyme-substrate complex is formed reversibly, while its dissociation into the 
product is considered irreversible. Further, the phosphate donor, which typically would 
be ATP, is assumed in constant concentration and not modeled explicitly. This reaction 
set-up is frequently used for studying signaling cascades; see e.g. [91 [Ft ] [T8 l [2? 1 [29 l I3T]. 

Assuming mass- action kinetics, the differential equations describing the dynamics of 
the system over time t are given by: 



dSl 
dt 
dS? 



(^i+^i)^% + W + ^ 1 



— = b°Y° 
dt 1 1 



dt 
dYl_ 

dt 
dY[_ 

dt 



= {b\+c\)Y l 1 -a\F l S}, 
= -{b° i +4)Y? + a < lS<lSU, 
= -{b\+c\)Yl+a\F l Sl 



[a i+1 S i+1 + a i Fi)S i , 



(2.2) 
(2.3) 
(2.4) 
(2.5) 
(2.6) 



for i = l...,n. Here a*,b*,c* are positive rate constants. For convenience we put 
for j = 0, n + 1, such that equation (2.2) also holds for Sq = E (i.e. 



a j 
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Figure 1: One-site PTM cascade of length n. The enzyme mechanism follows the classical 
Michaelis-Menten model. The reaction constants 8i,ji, A« are introduced: Si is the inverse 
of the Michaelis-Menten constant for Fi, Xi is the relative catalytic efficiency in layer i, 
and 7i is the catalytic efficiency of Fi divided by the dissociation constant of Sl_ x , c°. 



for i = 0). The steady states are found by setting the right hand side of the above 
equations to zero, which results in a system of 5n + 1 polynomial equations in 5n + 1 



variables with coefficients in R (for given rate constants). However, (2.4) and (2.6) result 
in the same equation and similarly other equations are redundant because of the following 

conservation laws: 



Fi = F i + Y?, E = E + Y 1 °, S i = S? + Sl + Y9 + Yl + Y? +1 , 



(2.7) 



for i = l,...,n, and Y^ +1 = 0. These can easily be verified by differentiation. The 
quantities E, Fi and Si are called the total amounts of enzymes and substrates, or just 
the total amounts. Note that E can be seen as the total amount of a 0-th layer with 
only modified substrate. 

From the form of the conservation laws we see that the steady state equations for Sf 



and Yl ((jOJ) and (j^), and that corresponding to E (p^, for i = 0), are redundant. 
In addition, the steady state equation for Y® +1 (2.5 
c i Y i > + b l Y i~- a i F i S i = °> which together with 



implies that (2.2) can be rewritten as 
2.4[) and (1231) reduces to c^YP-cjY^ = 0. 



Therefore, the system of equations (2.2), (2.4) and (2.5) is equivalent to 



0, (b)Yl 



at 



H + c t 



T F t Sj = 0, (c) Yi 



o 



6? + c' 



for i = l,...,n. Further, this system is equivalent to the system {(a) + cj(b), (6), 
(a) + 4(6)+c°( C )}: 



Y^- li F i Sj = 0, Yl 



5 l F l S} = 0, XiFiS} - S?S} 



i-1 
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for i = 1, . . . ,n, with the constants defined in Figure [TJ The constant 5i is the inverse 
of the Michaelis-Menten constant for Fi, Xi is the relative catalytic efficiency in layer i, 
and 7i is the catalytic efficiency of Fi divided by the dissociation constant of S'?_ 1 , c^. 



The system consisting of (2.7) and (2.8) is called SS (for steady state) and the so- 



lutions to it are the steady state concentrations for given total amounts. Let Con = 
{'jijSi, Aj| i = 1, . . . , n} U {Si, Fi, E\ i = 1 . . . , n} be the set of parameters of the system. 
Then SS is a system of polynomial equations with coefficients in IR(Con) and with vari- 
ables being the chemical species of the cascade. It consists of 2n quadratic equations (the 



last two equations in ( 2.8 ) for i = 1, . . . , n) and 3n + 1 linear equations (the conservation 



laws (2.7) and the first equation in (2.8) for i = 1, . . . , n). We are going to show that 
for fixed total amounts, the steady states can be found from the zeros of a polynomial 
in one variable. 

We note the following useful lemma: 

Lemma 2.9. At any steady state and for i = 1, . . . , n, 

(i) If E = 0, then E = 0. 

(U) If Fi = 0, or 1 + SiSj = 0, then Ft = 0. 

(Hi) If Sj = 0, then either E = 0, or Sj = for some j < i. 

(iv) If Sf = 0, then either Fi = 0, or Si = 0. 

If we allow i = 0, then statement (i) is included in (Hi). From here onwards we 
assume the total amounts are positive. 

By Lemma 2, Sj =0 and 1 + 5iSj =0 are not solutions to the system SS for any 
positive total amounts. Therefore, we can express the steady state concentrations of 
Fi, Yj 1 , Y®, and Sf in terms of those of the modified substrates: 

tp _ F i x _ SjF jSj Q _ nF jSj Q _ XjF jS} 

1 l + hsr Yi i + Yl i + 8 t s}> bi (i + SiSDsU [Z - W) 



The expression for Fi is found from the second equation in (2.8) and the conservation 



law of Fi. The other equalities are obtained by substitution of this expression of Fi in 



(2.8). We should remark here that the steady state values of YJ ,!^ 1 , and Fi depend 
only on the rate constants in the i-th layer, and on Fi and Sf. Additionally, they are 
expressed as rational functions of Sj. The steady state value of Sf, however, depends on 
the steady state values of the modified substrates in the i-th and (i — l)-th layers. 

2.3 Steady states as the zeros of a polynomial 

Let Xi = Sf for i = 0, 1, . . . , n, and x_\ = 1. By substitution of the expressions of K°, K 1 



derived from (2.8) into Fi,Si, the system SS is equivalent to the following system of 
equations: 

Y°- li F i Xi = Yl-hFiXi = i = l,...,n (2.11) 

XiFiXi- Sfxi^ = (l + 8iXi)Fi-Fi = i = l,...,n (2.12) 

Sf + Xi + (71 + 8i)FiXi + ji + iFi + iXi + x — Si = i = Q,...,n (2.13) 
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with the conventions So = E, F n+ \ = 0, and 5q = 70 = Aq = 0. Since the system admits 



no solution of the form x% = or l + 5iXi = (Lemma 2.9), we can eliminate the variables 



Fi and Sf in equation (2.13) using (2.12) to obtain an equivalent system given by (2.11), 



(2.12), and the equation 



AjFj(l + 5 i+ ix i+1 )xi + (1 + SiXi)(l + 5 i+ ix i+1 )xi-i(xi - Si)+ 

(7i + 5i)Fi(l + 5 i+ ix i+ i)xi-iXi + 7 i+ iF i+1 (l + 5iXi)xi-ix i+1 = 0, 

for i = 0, . . . , n. The left hand side of the equation can be written as a polynomial in 
Xi-i, Xi, Xi + \ with coefficients in M(Con), 

(ej-Xj+i + ef)xi + Xj_i(xj + i(efxi + efxj + ef) + efxf + ejxi + ef), 

for 7 = 0, . ..,n, with coefficients: el = Aj^+iFj, ef = AiF l; ef = <SA+i, ef = S l+ ieJ + 
ji+iSiFi+i, ef = -Si+iBi + ji+iFi+i, ef = S i: ej = 1 + {S l + ~f i )F l - 8 l S il and ef = -Si. For 
7 = 0, these reduce to e% = 5 l+ i, e§ = -8 i+1 E + j i+ iF i+ i, el = 1, and e§ = — E, with the 
remaining coefficients being equal to 0. 

Therefore, the initial system, SS, is equivalent to a system consisting of two sets of 
equations: a system, SS/, of linear equations in F^S^Y- and K with coefficients in 



M(ConU{xo, • • • , x n }), corresponding to equations (2.11) and (2.12), and a non-linear 
system, SS X = {qf = 0}, of n + 1 equations in n + 1 variables, X{, with coefficients in 
M(Con). 

This procedure exemplifies a general property of the steady states of signaling cas- 
cades. Signaling cascades are studied in full generality in [7], where we extend previous 
results of Thomson and Gunawardena [29]. We show that the steady states of a sig- 
naling cascade are determined by solving a system of M polynomial equations (with 
coefficients depending on the reaction rates and total amounts) in M variables. In our 
case M = n + 1, the variables are Sf, i = 0, . . . , n, and the system is SS X . 

Reciprocally, every solution to the system SS^ satisfying x% 7^ and 1 + 5{Xi 7^ for 
all 7, provides a solution to SS. The values of the other steady state concentrations are 
found from the linear system SS/. In other words, if we let 

n 

H(x , ...,x n ) = Y[ XiO- + 

then we have proved that the steady states of the system are in one-to-one correspondence 
with the solutions to the system SS X satisfying H / 0, that is, the solutions X{ lying 
outside the hypersurface Vh '■= {H = 0}. 

Remark 2.14. One can see that for rate constants and total amounts in general position, 
all solutions to the system SSa; lie outside Vh- Indeed, Xi = is a solution to the system 
if and only if h\ or h\ vanish, where h\ = (7i+2 + Si+2)Si — Aj+2(5j+i5j — 7^+1 Fj+i) and 
h\ = 7«+3 + <5j+3 — <5j+2Aj+3. The same conditions are obtained for solutions of the form 
1 + 

^3 x j ~ 0- I n this wa Y) only if the parameter set fulfils TJ i h\h\ — has the smaller 
system SSa; non-valid solutions. 
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The equations qf = may be used iteratively to eliminate variables. Indeed, note 
that qf can be expressed in the form 

qf (xi-i,Xi,x i+ i) = ai(xi,x i+ i) + a; i _ib i (a; i , 

whi ch is a linear polynomial in x\-\ with coefficients in M(ConU{xj, Xi+i}). By Lemma 
the solutions to SS satisfy 1 + <5jXj,Xj / 0. Hence aj(xj, Xj+i) = AjFjXj(l + 



2.9 



5i + iXi + i) / 0, and thus also bi(x{, Xj+i) 7^ 0. 

For i = n, we have g„(x n _i, x n ) = a n (x n ) + x n _ib n (x n ). This equation is used 
to eliminate the variable x n _i from = to obtain a new equivalent equation 

gf_ 1 (x n _ 2 ,x„) = a„_i(x n ) + x n - 2 K-i{xn) = 0, satisfying a n _i,6 n _i / for any so- 
lution to SS. We proceed iteratively in the same way for i = n — 2, . . . , 1 and obtain new 
equivalent equations, 

qf(xi-i,x n ) = ai(x n ) + Xi-ibi{x n ) = 0, (2.15) 

which are linear polynomials in Xj_i with coefficients in IR(ConU{x n }). Additionally, 
the coefficients are non-zero for any solution to SS. 

Finally, for i = 0, we find the equation q$ = (1 + <5iXi)(xo — E) + 71F1X1 = 0. 
Elimination of xo and x\ using <fi and q\, respectively, gives a polynomial in x n : 

(ToW) = (&2(x n ) - 5ia 2 (x„)) (ai(x n ) + £'&i(x n )) + 7i-Fi6i(x n )o 2 (x n ). (2.16) 

To sum up, we have shown that any solution to the system SS provides a solution to 
the system SS X = {qf = 0} together with SS;, and in particular that the steady state 
values of x n are roots of the polynomial q l Q{x n ). 

The reverse might not be true depending on the parameter values. Specifically, a 
root of the polynomial §q is a solution to SS if and only if cTj, bi, 5i-{di — bi 7^ 0. When 
this is the case, the values of xo, ■ ■ • ,x n _i can be found from the root x n and qf +1 as 
Xi = — aj + i(x n )/6j + i(x n ). The rest of the steady state concentrations are obtained by 



solving the linear system SS/ using (2.10). 



Note that equation (2.15) gives Xj_i = S}_ 1 as a rational function of x n = S^. These 



functions will be studied in detail in £2.7 



Theorem 2.17 (Polynomial). Consider a one-site PTM cascade with n layers. If posi- 
tive total amounts are given, then any steady state value of is a root of the polynomial 
Qoi^n)- Additionally, all other species concentrations are expressed as rational functions 
of Si 

Remark 2.18. The elimination of variables from the system SSa; could have been started 
from the equation for i = 1. This would have led to a polynomial in describing the 
steady states of the system. 

Remark 2.19. The relation Fj = F;/(l + 5 l S}) gives S} = (F l - F^/faFfi. Using 
this relation, the steady state values of S®,Sj,Y®, and Yf can be expressed as rational 
functions of F{ instead. It follows that the steady states can be described by means of a 
polynomial in F\ or F n . 
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We have proved that the steady states are the roots of a polynomial. Additionally, we 
have shown that the different variables are related through rational functions. Rational 
functions have quite well-described behaviors and this will enable us to derive theoretical 
results about the system. 

We are only interested in biologically relevant systems for which, in principle, it is 
possible to form all substrates and intermediate complexes. This suggests the following 
definition: A Biologically Meaningful Steady State (BMSS) is a steady state for which all 



total amounts are positive and all species concentrations are non-negative. Lemma 2.9 
ensures that for a BMSS all species concentrations are positive, i.e. non-zero. We call 
the set of BMSSs for D C M^ n+1 . 

It has been reported that multistationarity (more than one positive solution to the 
steady state equations for fixed total amounts) occurs in the MAPK cascade [20|, 123] . 
This cascade does not fit into our setting, because it has multiple site modification cycles 
in some layers. We will show in the following subsections that multistationarity cannot 
arise in our setting and that one-site PTM cascades have exactly one BMSS. According 



to Theorem 2.17 the value of determines the other concentrations at steady state. 
Therefore, multistationarity can only arise if the polynomial in has more than one 
zero for which all other concentrations are in D. We will see that there is exactly one 
valid root, which is the first positive root of the polynomial. 

2.4 Splitting the cascade into two parts 

We noted in the previous section that, provided positive total amounts are given, the 
steady state concentrations of a cascade of length n are the solutions to the system SS 



given by (2.7) and (2.8). 



Let us focus on the cascade obtained from the first i layers only. Its connection to the 
last n — i layers is through the intermediate complex Y® +1 , accounting for the conversion 
of Sf +1 to Sj +1 via the kinase Sf. If the steady state value of Y® +1 is known, then the 
steady state values of the chemical species in the first i layers satisfy the steady state 
equations of a cascade of length i with total amounts E, F%, . . . , Fi, Si, ... , S^-i, and 
Si-Y^^Sf + Sl + YP + Yl. 

Similarly for the cascade consisting of the last n — i layers. If Sj is known, then 
the steady state concentrations of the chemical species in layers i + 1, . . . , n, satisfy the 
steady state equations of a cascade of length n — i with total amounts i^i+i, . . . , F n and 
Si+l , • • • , S n . 

Therefore, the study of a cascade of length n can be reduced to the study of smaller 
cascades. This observation is key in what follows and will be used to derive most of the 
conclusions about the cascades using induction arguments. 

2.5 The system with only one layer, n — 1 

We first focus on understanding cascades of length n = 1, that is, a system consisting 



of only one cycle. According to (2.10), a steady state value of S 1 determines all other 
concentrations at steady state. 
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The system consists of the chemical species S , S , Y°, Y 1 , E and F. Using (2.10), 
the total amount S in (2.7) (with the last term now being zero) is 

„ AFS 1 . -fFS 1 5FS 1 

+ S +T^T + ^^T- (2-20 



(l + SS^E 1 + 5S 1 1 + 8S 1 ' 



Further, use (2.7) and (2.10) to express E in terms of E and S 1 , and conclude that S is 
a rational function of S 1 of the form 

5 = ^> = iW)' (2 ' 21) 

where A = 7 F - SE^qiiS 1 ) = 1 + 5S\ q 2 {S 1 ) =E — AS 1 , and piS 1 ) = XFqi(S 1 ) + 
qi(S 1 )q2(S 1 ) + ( , y+5)Fq2(S 1 ) is a polynomial of degree 2. The function </? can be extended 
continuously to S 1 = for which ^(S 1 ) = 0. 

Let r C 1 + be defined as T = R + if A < 0, and V = [0, E/A) if A > 0. 

Proposition 2.22. Let a one-site PTM cascade with n = 1 layer and positive total 
amounts be given. Then, 

(i) The system has a unique BMSS. 

(ii) There exists an increasing continuous function tp : M + — > T such that S 1 = t^(S) is 
the BMSS value of S 1 and with ip(0) = 0. 

( Hi ) Further, if) is the inverse of ip on V C ]R + . 

Statement (i) is well known and has been proven for instance in [H I30|. 132] . We note 
that ip' is always positive and thus ip is increasing in each of the intervals defined by 
the zeros of its denominator. Hence, it transpires that multiple steady states could only 
occur if A > 0. The previous proposition ensures that in this case the BMSS is located 
between zero and E/A and values of S 1 larger than E/A result in biologically non- valid 
concentrations. 

We now introduce another rational function relating E to Y°. If we isolate E from 
(2.20), we obtain E as a function of S 1 and this function does not depend on E. Sub- 
stituting S 1 = Y°/('-fF — 5Y°) (obtained from (2.10)) into this expression, we get the 
following relation at steady state, 

E = 9(y°) = Pl{ ^r v (2-23) 

where 

Pi(y) = Ay(£ - y), p 2 (y) = (5 + 7 )y 2 - 7 (l/«5 + F + £ + S)y + ^S, 



and £ = jF/5. 
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Proposition 2.24. The polynomial p2 has two positive roots a\ and a.2 with ot\ < £ < a<i 
and the function g is increasing and continuous in the three intervals defined by these 
two roots. The BMSS exists for values ofY° in (0,ai) only. Additionally, g(0) = 0, and 
g(Y°) tends to +00 as Y° tends to a\. 

It follows that the BMSS value of E is given as a continuous increasing function of 
the BMSS value of Y°. A choice ofE = E + Y° = g(Y°) + Y° imposes the BMSS value 
of Y°. 



2.6 One-site PTM cascades admit exactly one BMSS 

To proceed to higher values of n, we use an induction argument. Given a cascade of 
length n we consider separately the last layer C n of the cascade and the cascade C 
consisting of the first n — 1 layers. According to §2.4[ the link between the two cascades 
is the complex Y®. 

If we focus on the last layer C n , the kinase of this system is <S^_i- Recall that the 



function g n = g defined in (2.23) does not involve E. Thus, for any BMSS, the relation 



holds with constants j n , 5 n , X n , F n , and S n . By Proposition 2.24 g n is a continuous 
increasing function in [0, a%) with image K+. 

We focus now on the cascade C. As noted previously, C corresponds to a cascade of 
length n — 1 with total amounts E, F\, . . . , Si, ■ ■ ■ , S n -± — Y®. The key to proceed 

is that for a cascade of length n — 1, S^—i is an increasing continuous function ip n -i of 



5 n _i with V'n-i(O) = (this is proved by induction in Theorem 2.25). Therefore, since 



the total amount of substrate in the last layer of C is S n -\ — Y^, we have that 

SLi = f(Y°) :=Vn-l(5 n -i-F„°), 

with / a decreasing function in [0, SVi-i)- 

Thus, the split of the cascade provides two functions of Y® describing the BMSS 
value of S^—i- Both functions are continuous, g is increasing and tends to infinity as 
ai is approached, while ip n -i is decreasing. It follows that they intersect in exactly one 
point, which is the BMSS value of Y®. This is illustrated in Figure [2J The technical 
details of this argument are covered in the proof of Theorem |2.25[ 

Let A n = j n F n - SnSn-i = 5 n (£ n - S n -i) with £ n =j n F n /8 n such that Ai = A and 
£1 = £. Further, let T n = M + for A n < 0, and T n = [0, 5 n _i/A n ) for A n > 0, such that 

r 1 = r. 

Theorem 2.25 (One steady state). Let a one-site PTM cascade with n > 1 layers 
and positive total amounts be given. Then, the cascade has a unique steady state in D. 
Further, if S n is considered variable, then 

(i) Y® is a continuous increasing function of S n with image [0,/j, n ), where fi n = 
min(£ n , S„-i). 
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Determination of the BMSS 



°n-l 




as substrate 
(blue curve) 



&n-l — 9n(Xn) 

as kinase 
(red curve) 



S" 1 



F n 




BMSSs are located in the highlighted re- 
gion (before the first singularity of g n ). The 
two curves intersect in the BMSS value of 

(Xn > "Sn-l)- 



Figure 2: Induction argument to show the existence of exactly one BMSS 



(ii) S\ is an increasing continuous function ip n of S n with V'n(O) = 0. It is the inverse 
of a rational function <p n defined on . For n — 1, — (p and — ip. 

Remark 2.26. Our approach is not only useful for showing that there is exactly one 
BMSS. It also provides a constructive and iterative way of finding it. A priori we have 
no explicit analytical expression for vp n , but we do for its inverse (p n . 

Remark 2.27 (Effective computation). The rational function (p n and the set T n can be 
determined from the total amounts and the rate constants. Programs like Mathematica™ 
allow for the effective computation of an inverse function x = g(y) of y = f{x) using the 
command Interpolate. It proceeds from a list of pairs (y, g(y)) = (f(x),x). By generating 
many pairs (f(x),x), the estimated inverse function can be as accurate as desired. 

Remark 2.28 (Stability of the BMSS). We have shown that a one-site PTM cascade 
admits exactly one BMSS for any set of rate constants and any specified positive total 
amounts. This steady state is biologically attainable if it is (asymptotically) stable, that 
is, nearby trajectories are attracted to it. For a system of ordinary differential equations, 
a steady state is asymptotically stable if all eigenvalues of the Jacobian have negative 
real parts \34\ Thm. 1.1.1]. For n = 1, it is shown in pQ that the steady state is a global 
attractor. For 2 < n < 10, we randomly generated parameters ranging from 1CP 6 to 
10 4 . In all cases we have found that the stability criterion is satisfied and hence that the 
BMSS is stable. However, we have been unable to provide a mathematical proof. 

2.7 Rational functions relating the substrates 



In this subsection we are going to take a closer look at the rational functions relating the 
modified substrates derived from equation (2.15). We will provide an iterative expression 
for them that elucidates different biologically relevant properties. 
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For convenience, define the rational function relating Y® and Sf in (2.10) by 



Y? = 9 J{Sl) = ^f^. (2.29) 
Consider a cascade of length n. As we have seen, steady state concentrations in the last 



layer satisfy the relation Sj l _ 1 = g n (Y®). Combined with (2.29), this gives the relation 



ci — f ( qi\ — ^nF n Sn 

n - 1 ~ In - l{ J ~ dniSlO)' 

with 

di(x, y) = (Si - y) - x - F^Si + n)x + Si(Si - y)x - S^x 2 , 

1 < i < n. (The reason for introducing di as a function of two variables will be clear later.) 
In d n (x, 0), the sign of the leading coefficient divided by the sign of the independent term 
is negative, hence the polynomial has exactly one positive real root, say a n . For x = 0, 
we have d n (x, 0) > 0. Consequently, Sj l _ 1 is only meaningful for S [0, a n ) (is negative 
otherwise). Therefore, F n , S n and the rate constants of layer n set an upper bound 
to for any steady state, independently of the parameters in the previous layers of 
the cascade. In fact, the previous layers further constrain S^, e.g. S n -i establishes an 
upper bound to Sj l _ 1 and hence the maximal value of is further decreased. (Here and 
elsewhere, 'maximal value' is the smallest upper bound, i.e. the supremum, but it is not 
attained.) 

In general, consider the i-th. layer of the cascade. For every value of YR.i, the steady 
state values of the species in the first % layers are found by solving the steady state 
equations for the cascade consisting of layers 1 to i with the total amount of substrate 
in layer i being Si — Y®,y Therefore, applying the same argument as above, we have 

= 9i(S}, Yi+i) = , , cl J r j (2.30) 

with g n (Sn, Yn+i) = fn-i{S\), since Y®, x = 0. Arguing again as above, the polynomial 
di(Sf ,Y® +1 ) has exactly one positive real root in Sj, say Pi(Y® +1 ), for a fixed value of 
Y^e^Si). 

Consider < y% < y^. Writing yi = y\ + (y2 —yi), we find that 

di(p(yi),y2) = ~{V2 - yi) - Si(y 2 - yi)p{y\) < o 

and therefore p(y2) < p(yi)- In particular, if a« = Pi(0), then pi(Y® +1 ) < cti and Sj < a, L 
at any steady state. 

It follows that for every i, there is an upper bound, aj, for Sf at steady state, 
depending only on the rate constants and the total amounts, Fi and Si, in the layer. 
How far the steady state value is from the upper bound depends on the amount of 
sequestered substrate in Y® +1 . 
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Proposition 2.31 (Rational functions). For i = 0, ... ,n — 1, the BMSS value of Sj 
satisfies Sj = fi(Sj), where fi is an increasing rational function of Sj l defined on an 
interval [0,/%). It depends on Fj,Sj, j > i + Furthermore, < and < (3 n -i = 
a n for i < n — 1. In fact, the following relation holds 



A 



i+i ■ ■ ■ 



An-Fj+l • • • FnSn 



d i+1 (s} +v Y° +2 )---d n (si,oy 



where Sj 



M), j >i + l andYf = fJ{S l n ) 
recursively via (2.29) and (2.30). In addition, 



gj(fj(sl)), j>i + 2, 



are given 



(i) f3i is the first positive singularity of fi and a root of d- L+ i(x, Y® +2 ). In particular, 
the image of fi over [0, /%) is M+. 



(ii) If Si+i tends to +oo, then j3i tends to /3j + i (with j3 n 
(in) Sj < Qfj. 



+ooj. 



In Figure |3} we show the graphics of the rational functions, and the values and on 
for a cascade of length n = 3. 

The rational function for Sf_ 1 in terms of is obtained by considering S^-^ as the 
enzyme in the i-th layer. Therefore, we require layers 2, . . . ,n to compute the rational 
functions of S\, . . . , Sj^ in terms of S^, and layer 1 to compute that of E = 5q. Hence 
they are independent of E. When we require the rational functions to satisfy the con- 
servation law for E, the steady state value of S^ is determined and those of the other 



variables can subsequently be derived using (2.10) and Proposition 2.31 



From these results, one derives the following theorem. 

Theorem 2.32 (Downstream variability). Assume a one-site PTM cascade of length n 
is given with all total amounts but Si, for some i, fixed. Then, increasing Si causes the 
steady state values of Sh Yf , and Y 1 to increase for all j = i, 



,n. 



By Theorem 



2.25 ii), if S„ increases or decreases, then so does S„. Thus, if we 



consider the first j layers of the cascade, then an increase or decrease of the total amount 

On the other hand, by Theorem 



5 



Y® +1 will induce the same effect on Sj. 



2.32 



if 



all total amounts but Si are fixed, an increase in Sj causes the steady state value of YP 
to increase and therefore Si-i — Y® to decrease. It follows that Sf_ 1 and Y®_ t decrease 



(Theorem 2.32 ) . Proceeding with the same reasoning, we see that the effect is transmitted 
upstream in an alternating fashion. 

Theorem 2.33 (Upstream variability). Assume a one-site PTM cascade of length n is 
given with all total amounts but Si, for some i, fixed. Then, increasing Si causes the 
steady state values of Sj , j = 1, . . . , i, to increase for j = i — 2k, and to decrease for 
j = i-2k-l. 
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Rational functions 



We illustrate the rational functions and bounds in Proposition 2.31 for a cas- 
cade of length n = 3. Graphics (a-c) below depict the rational functions fi 
relating the steady state values of iS^, Si, E to S\. The first singularity, is 
highlighted in orange and delimits the region for the BMSS (in green). The a 
priori maximal value (the positive root of di(x,0) ) is also depicted. Note 
that a-i reduces the possible maximal value of S^. The highlighted region in 
graphic (d) corresponds to the inverse of the stimulus-response curve, see { 2 





a 2 = 3.77 






f /? 2 = 5.58 





(a) Si = / 2 (5 3 1 ) 



si 







i, = 2.47 








2 4 




r 


- 

A =4.27 







(b) ^ = /i(s!) 




So = 3.61 
(c) E = /o^ 1 ) 




cr, = 3.61 

(d)E = E + Y 1 ° = r(S 1 3 ) 



Parameters: 
Polynomials: 



diCic, o) 



= 1, Fi = 5, F 2 = 4, F 3 = 
■ 2a; - 0.5a; 2 , daO, 0) = 9 - 



3, Si = 8, S 2 = 9, S3 = 
■ 0.5a; - 0.5a; 2 , d 3 (x,0) - 



10 

= 10 - 



Figure 3: Rational functions for the modified substrate and stimulus-response. 



2.8 Stimulus-response curves as rational functions 



One-site PTM cascades have exactly one BMSS for any initial total amounts. Hence, for 
every value of E and all other total amounts fixed, there is a unique steady state value 
of 5*. The plot of against E is usually called the stimulus-response curve. 



By (2.29) and (2.30) we have 



E = E + Y 1 U = gi (Si,Y 2 ") + gY(Si), with 



n Y(ol) _ TiFiSj 



2.31). 



Note that fY = g\ ° fi * s a continuous increasing function of G [0, j3i) (Proposition 
As a function of Si, g\ is /o and hence by Proposition 



2.31 



it is a continuous 



increasing function of S\ £ [0,/3q). Since flo < Pi, we obtain the following theorem. 
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Theorem 2.34 (Stimulus-response curves). The stimulus-response curve = sr(E) is 
the inverse of a rational function 

E = r(Si) = = foiSi) + fYiSl), 

r 2Wn) 

for in [0, a n ), where a n = (3q is the first positive real zero of^iS^)- The polynomials 
r\ and r2 depend only on the parameters of the system and the specified total amounts of 
substrates and phosphatase (not the kinase). In addition, r is continuous and increasing. 

The plot in Figure [3^d) shows that r might have many singularities. However, the 
biologically meaningful region is given as the interval from zero to the first singularity. 

By construction, the polynomial q$ (S^) in (2.16) is exactly r2(S^)E — ri(S^) (up to 
multiplication by a real number). Therefore, we obtain the following corollary. 

Corollary 2.35. Let a one-site PTM cascade with n layers and positive total amounts 
be given. The BMSS value of is the first positive real root of the polynomial §q in S^, 
as given in Theorem 2.11. 

To finish the subsection we note the following. With fixed total amounts of substrates 
and phosphatase, all steady state concentrations are given as rational functions of S\, 
independently of E. This implies that the stimulus-response curve showing the response 
of any species in the cascade admits a rational algebraic parameterization in the form 

r c :[0,O -> R 2 + 

s H- {r(s),fc(s)), 

where C is a chemical species and fc(S^) is the corresponding rational function of 
for that species. The maximal values of Si, . . . , S^ l _ 1 are obtained as cjj = fi(o~ n ). 



3 Biological implications 

In section ^2] we have focused on the analytical description of one-site PTM cascades of 
arbitrary length n. The existence of exactly one biologically meaningful steady state has 
been shown. Additionally, our approach has provided explicit rational functions relating 
concentrations of substrates and enzymes at steady state. 

In this section we exemplify how our method can be used to provide qualitative insight 
into different biological aspects of signaling cascades. Generally, it is difficult to obtain 
experimental data and to estimate rate constants. Rate constants are typically only 
known in specific experimental contexts or to be within a certain range. It is therefore 
of importance to be able to derive conclusions that do not stringently rely on specific 
values of rate constants |13| 128] . We show that this is possible using our approach and 
that different cascade behaviors, for instance in response to noise or varied stimuli, can 
be studied using our analytical description. Further, it might be possible to design or 
guide experiments from the expected behavior of a system in order to uncover missing 
connections in a reaction network or the presence of feedback mechanisms |1 1| fT2| [28] . 

In what follows, stimulus refers to E while response refers to the steady state value 

of the modified substrate in the last layer. 
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Maximal response 

Reduction of the maximal value of re- 
sponse Sn for different cascade lengths 
n (for fixed parameters in the last 
layer) 

Length Max (<r n ) 
~1 5.58 

2 4.27 

3 3.62 

Parameters: 

a* = 6* = c* = 1, F 1 = 5,F 2 = 4, F 3 = 3 
Si = 8, S 2 = 9, S 3 = 10 



Response 




Stimulus 



Stimulus-response curves in semi-log scale (log(_E) vs 
Si) for different cascade lengths n — 1,2,3. 



Figure 4: Maximal values of the response S\ for different cascade lengths. 



3.1 Maximal response 

The maximal response of a cascade is the limiting steady state value of when the 
stimulus E is increased to infinity, or in realistic terms, it is the value of the response 
S^, when the stimulus is large. The total amount of substrate in the last layer sets an 
obvious upper bound to the maximal response. However, the maximal response can be 
far from the total amount of substrate. The first bound is imposed by the rate constants 
and the level of phosphatase in the last layer (by a n , cf. jj |2.7[ ). The bound is further 
restricted by the amount of substrate in the previous layer (since the substrate acts as 
kinase), and, as we will see below, the length of the cascade. 

We have shown that the maximal response attainable in a cascade of length n is given 
by the first positive zero of a polynomial, which corresponds to the denominator of the 



function r (Theorem 2.34). In Proposition 2.31 we saw that every layer of the cascade 



accounts for a reduction of the maximal response. In this sense, adding a new layer on 
top of the cascade lowers the possible maximal response. Only if the total amount in the 
new layer is very large, will the maximal response remain unchanged or essentially the 
same, cf. Proposition 2.31[ ii). 



It follows that each layer of the cascade reduces the maximal response by means 
of sequestration of substrate in the intermediate complexes. In Figure [4] this effect is 
illustrated for a cascade with n = 3 layers. The maximal response a n (= /So) of the last 
modified substrate is given in terms of the length of the cascade. In this system all rate 
constants are the same and the maximal responses are 5.58, 4.27, and 3.62, respectively, 
which are much lower than the upper bound set by the total amount (fixed to 10). 

3.2 Variation due to protein abundance 

Our setting is well suited to study how alterations in protein abundance in some layer 
(e.g. due to noise, or self-regulation) affects the system. The biological interest in this 
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Upstream and downstream variability 



a * = K = c * = !; F * = 3 . E = 3 > s « = 7 
Increase of S3 with different parameter condi- 
tions. It causes an increase of Sj , S3 , and S4 , 
and a decrease of Si . 




(a) Fs — 3. The substrate S3 goes 
to infinity and this causes Sg to ap- 
proach 04 = 3.27. 




(b) Fi = 10. In this case, S3 
does not tend to infinity and so 
S4 is kept low too. 



Figure 5: Illustration of Theorems 2.32 and 2.33 for n = 4. 



type of analytical inquiry is discussed in |18j under the name "slow regulation". It is 
pointed out that expression of phosphoproteins is altered during many biological pro- 
cesses. 

Variation in signaling protein abundance corresponds to variation in the total amount 
of substrate in some layer. For slow alterations, the cascade will readjust to a new steady 
state. Theorems 2.32| and |2.33| provide a clear picture of how modified and unmodified 
substrate concentrations of all layers are affected. 

Theorem 2.32 tells us that an increase in the total amount of substrate in some layer 
(say, i) is transmitted downstream at steady state as an increase in the concentrations of 
the modified substrates. Indeed, the effect is equivalent to increasing the stimulus in the 
smaller cascade consisting of the layers below the one undergoing variation. However, 
these layers have fixed total amounts and therefore the modified substrate in layer j > i 
is still bounded above by ay. As the total amount of substrate is increased indefinitely 
in the layer of modification, the modified substrate might also increase indefinitely or 
come to a halt at some specific value. In the latter case, the substrate accumulates in 
unmodified form. Whether one type of outcome or the other occurs, depends on the 
total amount of phosphatase in comparison with the total amount of substrate in the 
previous layer. 

Theorem |2.33 tells us that an increase in the total amount of substrate in layer 
i is transmitted upstream in an alternating way. In this case, the concentration of the 
modified substrate in a layer j above layer i depends on the variation in the total amount, 
Sj — Y^ +1 . Therefore, how large the change is, depends on the change in sequestered 
substrate itself Y® +1 . In particular, this suggests that the effect of variation is almost 
negligible in the layers far upstream of the varied one. This is illustrated in Figure [5] 

Regarding the concentration of unmodified substrate, then the effect is opposite to 
that of the modified substrate. By l&M = A^ / (7*^-1) > thus an increase in Sj 
causes Y® to increase and to decrease. Hence, Sf increases as well. However, for 
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j / i, the total amount Sj is constant and consequently, the concentrations S® change 
oppositely to that of Sj. 

3.3 Stimulus-response curves 

In our model, variation in the total amount of enzyme (stimulus) is a special case of vari- 
ation in the total amount of substrate in some layer. Stimulus-response curves have been 
the focus of many papers, e.g. [31 HU [181 EH EI] • A recurrent topic is (ultra)sensitivity of 
a system, or the capacity by which a signal transforms the output in a switch-like mode. 
An analysis of this requires a quantification of the system's switch behavior and different 
measures have been applied here jllj . 

Measures of sensitivity rely typically on the amount of enzyme (stimulus) required 
to achieve 90% of the maximal steady state value of the modified substrate (response) 
compared to that required to achieve 10%. Recall that we define o~ n to be the maximal 
value of the response For M in [0, 1] , let Em be the amount of E that corresponds to a 
steady state value of S\ = Ma n . For example, 90% of the maximal response corresponds 



to a steady state value of 0.9<r n . According to Theorem 2.34, the total amount of enzyme 
required is determined as E'o.g = r(0.9a n ). 

Therefore, our results provide a way to analyze and compute the sensitive or switch- 
like character of a signaling cascade. In particular, the response coefficient (also called 
cooperativity index) is given as R = Eq. 9/-E0.1 [H], the switch value as E'o.g — -E0.1 [TT] 
and the Hill coefficient as n# = log(81)/log(£'o.9/-E'o.i) [15] - 

Other measures of sensitivity, probably more suited for an analytical study, escape 
from our control because the analytical form of the stimulus-response function is required. 
In our case, the stimulus is given in terms of the rational function r evaluated in the 



response (Theorem 2.34), hence the stimulus-response curve is the inverse of r. We do 
not have an expression for the inverse and can only evaluate this through tabulation 
of points (Ma n , Em) = (Ma n ,r(Ma n )). One alternative measure is the control curve 
given as (x/ f){df /dx), where / is the stimulus-response function of the variable x |llj . 
For / = r , a plot of the control curve could be created similarly to that of r _1 by 
tabulation and using that the derivative is one over the derivative of r. 

Further, we have studied the shift in stimulus-response curves for modified substrates 
in different layers, as suggested in [3T| Fig. 6]. For that, we note that U{ = fi(cr n ) is the 



maxi mal value of the modified substrate in the i-th layer (Proposition 2.31 and Theorem 



2.34). Then, if M G [0, 1], we consider E M to be the amount of E that corresponds to a 
steady state value of Sj = Mai, that is, M times the maximal response. We have that 
E\.j = r(M fi(a n )) and for any set of parameters 

E 1 m<E 1 m (3.1) 

(proved in Appendix). That is, the stimulus-response curves for the different modified 
substrates are shifted from right to left with increasing index of the layer. Consequently, 
in order to achieve maximal response, 5„ requires the least kinase while S\ requires the 
most. This is illustrated in Figure [6} 
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Stimulus-response curve 

Shift of the stimulus-response curves 
of the different modified substrates Sj, 
i = 1,2,3, for n = 3. Substrate val- 
ues are normalized to allow comparison. 
Graphic in semi-log scale. 

Parameters: 

a* = ft* = c* = 1, Fi = 5,F 2 = 4, F 3 = 3 
Si = 8,S 2 = 9, S 3 = 10 



% Response 




Stimulus 



Figure 6: Stimulus-response curves in semi-log scale. 



3.4 Enzyme competition 

The one-site modification cycle considered here is driven by two opposing processes, 
modification and demodification, catalyzed by a kinase and a phosphatase, respectively. 
Both enzymes compete for the substrate and one expects that if the kinase is in excess 
over the phosphatase, then the system shifts in favor of the modified substrate. The 
results in section fallow for a precise statement of this fact, as well as its extension to 
the entire cascade. 



It follows from Proposition 2.22, that the steady state value of the modified substrate 
is an increasing function of the total amount of substrate. The modified substrate takes 
values in 1+ for A < 0, and [0,£/A) for A > 0, where A = ~/F-5E = ^((c 1 /c°)F -E). 
Hence A < only if T = c l F — c°E < 0. This implies that as the total amount of 
substrate is increased, the steady state concentration of the modified substrate increases 
indefinitely or approaches the asymptotic value E/A, depending on the initial amount of 



enzymes and the dissociation constants. Additionally, from (2.10) and the conservation 
law for E, we have that 5° = XFS 1 /(E — AS 1 ), which is an increasing function of S . We 
conclude that when S increases, S° tends to —XF/A > if T < and to +oo otherwise. 
Only when T = 0, both the unmodified and modified substrate concentrations become 
large at the same time. 

For cascades of length n > 1, the same competition is observed when the total 
substrate concentrations 5i,...,5 n _i are fixed and S n is increased. Theorem 2.25 
shows that Y® tends to min((c„/c°)F ri , S n -i) as S n becomes large. In addition, 



is a function of y„°, f(Y°) with /(£„_! ) = 0, and hence by ((^ and flgTO) , 5° 



ii 



(A„y n )/( 7n /(y„ )) and 5i = Y°/(j n F n - 5 n Y®). Therefore, the sign of T n 
c^Sn—i determines the limit value of Y®, and thus the asymptotic values of the modified 
and unmodified substrates. Indeed, if T < then S 1 tends to infinity, while 5„ tends to 
the asymptotic value; when T > the opposite behavior is observed. Only when T = 
both substrates become large at the same time. 
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4 Concluding remarks 

In this paper we have studied signaling cascades consisting of n identical layers of mod- 
ification. We have shown that it is possible to draw biological relevant conclusions from 
analytical aspects of rational functions describing the system. The complexity of the 
system is reduced by splitting the cascade into smaller cascades. 

In particular, we have shown that signaling cascades with one PTM in each layer 
cannot exhibit multistability for any rate constants or total initial amounts of substrates 
and enzymes. It is well known that more complex signaling cascades, like the MAPK 
cascade, can exhibit bistability, e.g. [HI 120^ [23] . Since multistability in PTM only arises 
when there are multiple modification sites, we conclude that multistability in a signaling 
cascade (without feedback) must be a consequence of the presence of multiple steady 
states in some of its layers, and cannot be an effect of the cascade itself. 

Further, we have shown that stimulus-response curves can be obtained as inverse 
functions of certain rational functions which are explicitly given. Stimulus-response 
curves are important as they provide theoretical interpretation of the behavior of the 
system. Also, points on a stimulus-response curve can be determined experimentally 
and used to draw inference on the kinetic parameters of the system, e.g. |21j . However, 
the stimulus-response curve cannot be given in the form of a closed analytical expression. 
To find the inverse of the corresponding rational function we need to determine the roots 
of a high-degree polynomial. For n = 3, the polynomial describing the steady states of 
the system has already degree 17 rendering exact analysis difficult. 

Finally, we have studied sequestration. Variation in the levels of the total amounts of 
substrates influence the steady state of the system. Downstream, variation is transmitted 
positively, while upstream, variation is transmitted in an alternating fashion. As a 
consequence, the modified substrate in the last layer is always positively influenced by 
variations in any layer. 

The framework presented here is quite flexible and allows us to address questions 
about the system theoretically without restoring to simulation. The way we use the 
modularity of the cascade to derive properties of the cascade and to determine the 
number of steady states suggests that our results could be extended to signaling cascades 
with layers consisting of other reaction systems, as long as we have some mathematical 
knowledge about the variation of the species with respect to the total amounts. Such a 
study should aim to relate the number of steady states of the cascade to the number of 
steady states arising from each isolated layer. As a consequence, a better understanding 
of the emergence of multistability in signaling cascades may be provided. In particular, 
we believe that the ideas developed in the paper could be useful for studying more 
complex cascades, like the MAPK cascade, and for example help to better elucidate the 
parameter regions of bistability [22] . 
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A Proofs 



Lemma{2 

then Y"] = (since E 



We use <\1^ and (|2 

Sq by convention) and hence E 



then = 0. If 1 + 5iS} = 0, then Y> 
then y.° = Yl = 0. Also either (a) S° 



repeatedly without further reference, (i) If E 

E + Y? = 0. (ii) If Fi 



-F; and hence Fj 



0. 



0, 

0, 

= 0. (iii) If S} = 0, 
If Sf = (a), then 
0, while for i < n, 
i _ 1 = (b), then repeat 

the argument until either Sq = E = (hence E = 0, according to (i) ) or S® = for some 
j < i. In the latter case it follows that Sj = 0, using (a), (iv) If S® = 0, then Y® = 0, 
and either Fi = or Sf = 0. The result now follows from (iii) (a). □ □ 



yO 



Sf + S} + YP + Yl + YP +1 



Y° 



or (b) SU = 
For i = n, Y°, 



(ji+l j Aj+i) S\ S® +1 = 0; hence in both cases S{ = 0. If S}_ 



Proposition 2.22. The function ip is well defined for S 1 = with value S = 0, and it is 
continuous outside the roots of the denominator. The only possibly positive root of the 
denominator is r/ = E/A, provided A^0. 

The derivative of ip with respect to S 1 is ip' = 1 + N/[qi(S 1 ) 2 q2(S 1 ) 2 ], where N = 
XFql^-q'zS^+eqliqi-q^S 1 ) ande = 7 + «5F. Since q 2 - q^S 1 = E and qi-qiS 1 = 1, 
then N > 0. Therefore, ip' > for all values of S" 1 for which tp is defined, and hence it is 
an increasing function. 

2^91 Since E > 0, we have E ^ at 



In the following we rely on (2.10) and Lemma 



steady state. Therefore in order for E > 0, we need 



E = E-Y U 



E 



jFS 1 
1 + 5S 1 



> E > AS 



(A.l) 



For values of S 1 for which (A.l) holds, we have ^("S" 1 ) > and qiiS 1 ) > 0, hence 



piS 1 ) > 0, and ^(S 1 ) = if and only if S 1 = 0. Therefore, S = ^(5 X ) > if ((All 



holds. Define the remaining quantities (S°, Y°, Y 1 ^) through (|2.10[). These quantities 



are all positive. Therefore, a positive value of S 1 provides a steady state in D if and only 
if (2.21) and (A.l) are satisfied. The question is whether for every value of S there is a 



unique value of S satisfying equations (2.21) and (A.l). 

Let Ip be the restriction of ip to the set T of values S 1 
two possible scenarios: 



> satisfying (A.l). There are 



(i) If A < 0, then (A.l) is satisfied for all values of S and so T 



L|_. In this case, 

there is no positive root of q2(S 1 ) and hence Ip is a continuous increasing function 
in r. 



If A > 0, then (A.l) is satisfied if < S < rj = E/A. Since the only positive 



root of the denominator of tp is rj, we have tp is a continuous increasing function in 

r = [o, v ). 



In both cases, the image of tp is the set of positive real numbers M+. By Lemma 2.1 
the inverse function, ■0, is also a continuous increasing function of S in M + with image 
set r. It follows that there is a unique steady state of the system in D for positive total 
amounts E, F, and S. □ □ 
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Proposition 2.24 It is straightforward to check that the derivative of g = p\jp2 as a 
function of Y° is always positive and hence g is increasing, whenever it is defined. Recall 
that Pl (y) = Ay(£ - y), p 2 (y) = {5 + 7 )y 2 - j(l/6 + F + £ + S)y + 7 £S, with f = jF/5. 
The signs of the coefficients of P2(y) indicate that pziy) has either two positive roots 
or two non-real conjugate complex roots. For y = 0, P2(0) > 0, while for y = £, 
f>2(£) = — 7£/<5 < 0. Therefore, there is at least one positive root a± £ (0,£). When 
y tends to infinity, P2(y) becomes positive, and thus the other real root a 2 satisfies 
a2 > £ > a i- It follows that P2(y) is positive in [0, ai) and (a2,+oo) and negative in 
(cti)«2)- Note that pi(y) has constant positive sign in (0, £) and takes negative values 
in (£, +oo). Hence, E is positive only if Y° £ (0, a%) U (£, 02). 

If Y° > £, then Y 1 > F meaning that F < 0. Therefore, any steady state in D must 
satisfy Y° G (0, £). It is left to the reader to check that the derived steady state values 
for Y 1 , S 1 , S°, F are in D for Y° G (0, ai). 

In the interval [0, a±), g is a positive continuous increasing function of Y , which is 
zero when Y° = and tends to infinity when the root a\ is approached. □ 

□ 

Theorem\EM We are going to prove the theorem by induction on n. For n = 1, (ii) 



2.22 



follows directly from Proposition 2.22 ; (i) follows by composition of ipi = %j) in Proposition 
with the expression for in (2.10). Next, assume that the theorem holds for n — 1 



for some n > 1. 

Consider the splitting of the cascade of length n into the last layer C n and the first 
Ti—l layers C. By induction hypothesis, for every value of Y® £ [0,S n -i), there is 
exactly one steady state of C. In addition, is given by a continuous increasing 

function ij) n -\ of the total amount of substrate in the last layer. Hence = f(Y®) 

with / a decreasing function of Y® in the interval [0, S n —i). Further, / can be extended 
continuously to S n -\ with /(<S n -i) = 0, since ^n-i(O) = 0. Note that / is independent 
of S n . 

On the other hand, we noticed in the main text that S^_i is given by a continuous 
increasing function g n of Y®, for Y® in the valid interval [0, a\). The function g n is zero 
when Y® = and tends to infinity when the root a\ is approached. 

Therefore, there are two continuous functions of Y®, f(Y®) and g(Y®), describing the 
steady state value of S^-i- Both of them take positive values, but f(Y®) is decreasing, 
while g(Y®) is increasing in [0, a±). In addition, /(SVi-i) = and /(0) > 0. Since the 
image set of g on [0, ai) is M + , there exists a unique value of Y® for which f(Y®) = 
g(Y®), and therefore there is a unique BMSS. Note that the intersection point satisfies 
< y n ° < min(£ n ,S n _i), f n = ~fnFn/5 n - 

What remains to prove is that the statements (i) and (ii) also are satisfied for n. (i) 
The steady state is the unique value for which f(y) = Pi{y)/p2{y)- Note that p 2 {y) is a 
linear polynomial in S = S n with coefficients in Therefore, P2{y) = a(y) + b(y)S 

with b(y) = ^ n {lnF n — S n y)/5 n > 0. The steady state solution is then given by 

g ( \ pijy) - a(y)f(y) 

s=My) = mm ■ 
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This function is continuous in the interval [0, /i n ), where /i n = min(£ n , S n -i). In addition, 
its derivative is positive. Indeed, 

= Pibf + f(ab' - a'b) + f(p[b - pib')) . 

We have that /' < and p 1: b > 0. Also, ab' - a'b = F^/5 2 + -y(5 + -y)(Fj/6 - y) 2 > 
(suppressing the subindices) and p\b — p\b' = 7A(£ — y) 2 > 0. Consequently, ip' n > 
and ip n is increasing. Therefore, by Lemma |2.1[ the intersection point is given by a 



continuous increasing function of S. It proves (i). (ii) From (2.10), 



1 o Y n 

s n = x(Y n ) = _ ™ , 

fn- 1 n l -'n 1 n 

which is an increasing function of Y®; hence by composition of % with ip^ , S n i s an 
increasing function of S n . The inverse of this function is p n {lnF n S n / (1 + 5 n S^)), which 
is a rational function of S n , as desired. It is defined on r n = [ 0, xiHn)), which is +oo if 



p n = £ n , and S n -i/ A n otherwise. The proof of Theorem 2.25 is completed. □ □ 



Proposition 2.31, The proof is by decreasing induction starting from i = n — 1. The 



discussion above Proposition 2.31 shows that there is a rational function S n _ t = f n -i(S^) 
satisfying (i) with /3 n _i = a n . To prove that the image is K+, note that the derivative of 
fn-i with respect to S n is always positive, and therefore S n _ 1 is a continuous increasing 
function of S n (for S n , F n fixed), which tends to infinity as S n tends to a n , and vanishes 
if S n = 0. (Note that f n -\ only depends on F n and S n: so when S n (at steady state) 
tends to a n , it implicitely implies that S n -i tends to infinity.) (ii) follows from the 
expression of d n (x,0), and (iii) is likewise fulfilled. 

Now assume the proposition is true for all j satisfying i < j < n — 1 for some i. Con- 
sider i — 1. By induction hypothesis, Sj = fj(S n ) are increasing rational functions for all 
j > i, and so, by (2.10), Yj* +1 = fj+i(S n ) is an increasing rational function of S n . Hence 
by ( |2.30[ ), is a rational function of S n given by = gi(fi(S^), fl +l (S^)) . If 

i + 1 = n, then f n is the identity function. Note that f n -i depends on S n and F n , but 
not on any other total amounts. Similarly, g$ depends on Si and F{. By induction /j_x 
only depends on Sj and Fj, j > i. 

By induction hypothesis, Sf, Sf +1 and YR t are continuous increasing functions of 
S^ defined on the interval [0, /%). In addition, as S n tends to Sj tends to infinity, 
Sf +1 tends to the value rj = /j + i(/3j) < +oo (because < /3«+i) and so Y® +1 tends to 

V Y = &i(A) <+<*>■ 

The function pi{Y® +l ) is a continuous decreasing function of Y® +1 , ranging from m to 

0, defined as the positive root of the polynomial di(x, Y^_ x ). Let ~p i = p^ofX v Since Y® +1 

is increasing as function of S n , then Pi(S n ) is decreasing, because Pi(Y® +1 ) is. In order for 

Sl_ x to be well defined, it must be that Sj = fi(S n ) < PiiS^) (it is negative otherwise). 

Therefore, the possible steady state values of S n must satisfy fi(S n ) — Pi(S n ) < 0. The 

function /j — p i is a continuous increasing function that takes the value — at S n = 

(Y® +1 = 0). On the other hand, when S n tends to /3j, the function /j — p i tends to +oo. 
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Therefore, there exists a unique value of S^, ft-i < ft, for which f i — p i = 0. This value 
is the upper limit for Sj l . Thus, a necessary condition for valid steady states is that 
Sj l G [0,ft_i). Note that /i(ft-i) = pj(ft_i) is the Sj-v&lue corresponding to the first 
positive root of fKi(S^)) . 



It follows from (2.30) that Sj_ 1 is an increasing rational function of Sj 6 [0,ft_i), 
because the denominator is non-zero for values < ft-i- It proves the first part of 
the statement. Since for positive values of S*„ smaller than ft_i the denominator is non- 
zero, we conclude that ft_i corresponds to the first positive zero of the denominator. 
In addition, when Sj l tends to ft-i, the denominator of (2.30) tends to zero, while the 



numerator tends to some positive real number (since ft_i < ft) and hence S , /„ 1 tends to 
infinity showing (i). 

To see (ii), note that fi does not depend on Si, while p { does. If a is a fixed value of 
S^, and Si < Si, then we have Pi(a, Si) < J>i(a, 5J. Therefore, the intersection point of 
the two curves increases if Si does. In addition, pi tends to infinity as Si does (consider 
the expression for the positive root), and therefore intersection points as close to ft as 
desired can be obtained. □ □ 



Theorem ] 2. 32\ From (2.10) and Proposition 2.31 we have that Sj and Yj\ 1 are continu- 
ous increasing functions of Sj. Hence, by the inverse function theorem, Sj is a continuous 
increasing function of Y® +1 . Call this function ^j. It does not depend on Sj. 



On the other hand, the proof of Theorem 2.25 demonstrates that for every fixed total 
amount Si, Sj is a continuous decreasing function of Y® +1 , $i(Si,Y® +1 ). It follows from 
the fact that the total amount in the i-th layer, Si — Y^ +1 , decreases with Y® +1 . Therefore, 
the steady state value of Sj is the intersection of these two curves, \Pj(— ), $i(Si, — ). If 
now Si is increased, the first curve does not change, while the second one does. For every 
fixed value v of Yj%, we have that $i(Si, v ) < 3>i(«5j-, v) if Si < S^ Therefore, the steady 
state value of Sj increases as Si does, because we intersect with a growing function. 

The quantities Sj,j>i, are continuous increasing rational functions of S^. Further, 
for j > i, the rational functions do not involve Si- Hence, S^, . . . , Sj +1 will also be 



growing with Si. It follows from ( |2.10[ ) that Yf and Yj also increase, while Sj decreases. 

□ □ 

Theorem \2. 331 We proved that an increase of Si causes Sj to increase and S?_ 1 to de- 
crease. Using an induction argument, assume that the statement is true for layers with 
indices smaller than some value k. Then, if k = i — 2j, by induction hypothesis Sj_ 1 
decreases if Si increases. Therefore, by (2.10), Y®_ 1 decreases and hence Sk — Y^_ 1 
increases. It follows that Sj increases, as desired. If k = i — 2j — 1, the argument is 
analogous. □ □ 



Inequality (3.1). Since S^ is a continuous increasing function of Sj, we have that E = 
r % (Sj), where r l = r o /r 1 , is given by a continuous increasing function. Additionally, 
Sj is a continuous increasing function of Sj +1 given by g i+ i = fi° on [0, cri+i). We 



obtain r l+1 



r o 



g i+ i- Let E M = r % (Mo~i) be the value of E corresponding to Mai, the 
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maximal value of Sj. Because <7j = g i+1 (o~i+i), we have E M = r % (M g i+1 (ai + \)) . On the 

other hand, E^ 1 = r*+ l (Ma i+l ) = r^g^Mai+x)). 

Because r l is an increasing function, the proposition follows if g i+1 (Mai + i) < Mg i+1 (o-.i + i). 
Let Y < Y' be the values of Y® +2 corresponding to M<7j + i and cxj+i, respectively (zero if 
i = n-l). Since d i+ i(Ma i+ x, Y) > d i+1 (a i+1 ,Y'), then by (|2.30|), 



and the inequality is proved 



Ai+i-Fj+iMcTi+i \ i+ iF i+ i(Ji + i , 
1 — ni7 vn < M w — ^ = M 5i+i(^i+i, 



□ 



□ 
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